LiV204: frustration induced heavy fermion metal 
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We propose a two-stage spin-quenching scenario for the unusual heavy fermion state realized 
in the mixed valent metal LiV2 04. In this theory, local valence fluctuations are responsible for 
the formation of partially quenched, spin i moments below room temperature Frustration of 
the intersite spin couplings then drives the system to produce the heavy Fermi liquid seen at low 
temperatures. |Q] The anomalous resistivity and the sign change of the Hall constant can be under- 
stood naturally within this model, which also predicts a unique symmetry for the heavy quasiparticle 
bands that that may be observed in de Haas- van Alphen experiments. 

PACS numbers: 71.27.-|-a,75.20.Hr 



Magnetic frustration has attracted increasing attention 
in recent years as a new tool to suppress antiferromag- 
netism and generate new types of electronic behavior. 
One exciting possibility is the development of frustration 
induced heavy electron behavior in d-electron systems 
where the Kondo temperature is generally too small to 
overcome magnetic order without frustration. The recog- 
nition of heavy fermion behavior in LiV2 04[|| was an 
important first step in this direction. In LiV204, the 
frustrated lattice remains undistorted to the lowest tem- 
peratures measured, Q The heavy electron state which 
then develops is unusual in many respects: it displays a 
monotonically increasing resistivity [Q, a pressure-driven 
metal-insulator transition^ and field independent heat 
capacity up to 30TQ. 

LiV204 forms a spinel structure. The magnetic vana- 
dium atoms are homogeneously mixed valent forming a 
lattice of undistorted corner-shared tetrahedra. The spin 
at each site fluctuates between S = 1 (3d^) and S = ^ 
{3d}) with a formal valence of 3.5. Antiferromagnetic 
couplings between these spins Q imply that the lattice 
is a highly frustrated metal. There is a strong on-site 
Hund's coupling between spins of order 1 eV|^, which is 
manifested in the photoemission[^ and the high temper- 
ature magnetic susceptibility [Q. At room temperature 
the magnetic suceptibility and specific heat indicate the 
presence of spin i local moments at every site.[^ ||]. Be- 
low 4K, the moments quench into a heavy Fermi liquid. |^ 

There have been several attempts to explain this 
physics. One class of model sees the frustration as 
secondary, the essential physics being due to mixed 
valence JTot, exhaustion|^, ID chains |Tl|| or a two-band 
Hubbard model Another argues that frustration is 
essential, leading to a tetrahedron rule|]l3j, a frustrated 
Kondo lattice jl4j, antiferromagnetically coupled spins 
Hund's coupled to conduction electrons or quantum 
criticality|l^. Here we present a simple two-band model 
which we believe captures the essential physics: a two 
stage spin quenching that could not occur in the absence 
of frustration; and show how our results can be extended 
to the pyrochlore structure. 



The high temperature physics of LiV2 04 can be un- 
derstood as an S=l Mott insulator doped with holes on 
a pyrochlore lattice. The strong Hund's coupling leads us 
to consider only spin-1 and spin-i sites coupled antifer- 
romagnetically over the lattice via direct exchange |17|, 



H = -tJ2iXoa{^)Xao{J) + h.c.) + J2Js^ 

(i,3) (jj) 
,cr >< l,cr'| = \^h„ 



(1) 



where Xq„ — \\,<y >< 1, cr'| = x^^ct is a Hubbard opera- 
tor allowing the charged holes to move in the spin back- 
ground via a hopping matrix element t, subject to the 
constraint bl^b^ + x'^X = '^^ ^ is an antiferromagentic cou- 
pling between nearest neighbor spins. The local physics 
of this model is that of a mixed valent impurity under- 
going charge fluctuations d? ^ d} + via hybridization 
with the surrounding bath of holes, as might be formal- 
ized in a dynamical mean-fleld description of the above 
model. The effective impurity Hamiltonian replacing the 
first term in H would then be 



(2) 



As the frustration of the pyrochlore lattice prevents 
magnetic ordering from occuring, this term dominates 
the high temperature physics. From strong coupling 
arguments fl^, we know that a Fermi sea of electrons 
forms below Ti, quenching the local moment character 
on each site to S=^. This provides a natural explanation 
for the behavior seen below room temperature. On the 
lattice, this fiuid can be visualized as a smooth Fermi sea, 
from which the residual spin ^ moments protrude like 
spines on the hide of a porcupine. (Fig. 1.) The original 
antiferromagnetic interaction between nearest neighbor 
sites now couples the electron sea to the residual spin ^ 
moments. The properties of this "porcupine" Fermi sea 
are the topic of the remainder of this paper. The question 
is - what is the nature of the intrasite screening process 
that removes the free spins (spines) at low temperatures? 

At low temperatures, J becomes an effective Kondo 
coupling between the localized spins and the newly freed 
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FIG. 1: (a) High temperature phase with large moments, (b) 
Intermediate temperature phase with underscreened S = 1/2 
moments embedded in a Fermi sea hke the spines of a 'Porcu- 
pine' (left: corresponding j filled band interacting with a lo- 
calized level) (c) Low temperature phase in which the Heisen- 
berg interaction hybridizes these bands leaving a heavy Fermi 
liquid of holes, (left: corresponding band picture) 

electrons in the Fermi sea. The locally symmetric spin 
state required to satisfy the strong Hund's coupling 
means that spatial overlap of the realized Wannier states 
is forbidden as the heavy quasiparticle forms at low tem- 
peratures. Thus, even after the complete break-down of 
any orbital picture, as the mixed valent spin quenching 
of the upper level occurs in the totally symmetric chan- 
nel the second stage spin quenching can occur only in an 
orthogonal channel. This is the analogue of the atomic 
Hund's selection rule for itinerant states. Thus the con- 
ditions 



no overlap 

no phase transition (equivalent sites) 
no time reversal symmetry breaking 



arise, uniquely specifying the second channel (Fig. 2). 
The unusual symmetry of the composite quasi-particles 
resulting from this second spin-quenching should be ob- 
servable by de Haas - van Alphen experiments. 

We now treat the second stage at the mean field level 
keeping in mind that for LiV204, we really want to have 
a half-filled localized d-band interacting with a quarter- 
filled dclocalized d-band. Estimates of J « d^w'^ ~ 
400is:§, Tc2 « 4 - bm and t « 2eV^ (band theory) 
have been made, and we expect that Td (x t and Tc2 oc 

te J as these are the two energy scales in the problem. 
It is convenient to recast the Hubbard operators into a 
slave boson language. 



<y> 



(3) 



where Xq^ = b^Ca- is a Hubbard operator allowing the 



holes to move in the upper e^ symmetry channel, s 



4 4) 



is the total spin, and we need to en- 



force the constraints 



nb = 1, = 1 (a localized 



spin-i level of aig symmetry). The Hund's coupling Jh 
is to be taken to infinity to enforce the constraint S*^ = 2 
at doubly occupied sites. 

We now demonstrate how a mean field treatment of 
Eq.3 leads us to the effective Hamiltonian: 



Ck<j 



Ai(|6|2-1)-A2 



where 



-2i|6|2$ 



(1) 
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cl 
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(5) 



Clfc 



CNk) demark the 



number of bands arising from the magnetic unit cell of 
the lattice under consideration. Ai and A2 are Lagrange 
multipliers fixing the occupancies of the electron and spin 
liquid. is an N dimensional matrix representing the 
Wannier functions, where N is the number of atoms per 
unit cell and F = 1, .., M specifies the number of orthog- 
onal channels required to decouple the Heisenberg inter- 
action between one site and its M nearest neighbors. 
The mean field kinetic energy term of H is 



Hkin = 



<ij> 



ticibUb^c,), = -2t\b\^J2^l^k (6) 



where the bosons take an expectation value, and the 
fourier transform happens to give the totally sym- 
metric Wannier state. Magnetic frustration enables us to 
neglect the normally dominant interaction between the 
residual spin i moments, however, we must still con- 
sider the Kondo coupling between localized and dclocal- 
ized bands in the expansion of the Jsi ■ Sj term of Eq.3: 



H, 



<ij> 



where the operators of the dclocalized level have been 
expanded in terms of orthogonal linear combinations as 
-0|p — (f>*^{j — i)cj and (j)*^ is a form factor reflecting 
the phases at different nearest neighbor sites of i, reflect- 
ing the possible arrangements of c^ on the lattice. This 
becomes: 



E 



A»rA,r 



(8) 



where we have performed a Hubbard-Stratonovich gauge 
transformation on the four-fermion term before fourier 
transforming to obtain the off-diagonal terms in Eq.(5). 
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To treat the term of interest from a Hund's coupling: 
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H 



H 
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(cida^dlcp) ^ — + + (9) 



a Hubbard-Stratonovich gauge field has been introduced 
to decouple the four-fermion term|l^. Before the decou- 
pling, Hh\S = 1 OIS* = 1 > and Hh\S = >= 
^IS* = >, so in the limit Jh oo, the 15 = > state 
will be forbidden. To write Eq.(4) we have taken the limit 
Jh oo. Although this appears to be non-Hermitian, 
the mean field values of ^, ^ (independent variables) will 
also be imaginary so that physical quantities are real. 
The Free energy can then be expressed as 



F = -r^Tr[/n(l + 

+ Ai(|fe|2-1)-A2 



ArAr 



(10) 



where is the Hamiltonian expressed as a matrix as 
above. Since f and A both couple to Cl^Va, changing 
one affects the other, such that JA + g^f^ — 

0. Integration over such fluctuations yields -jr ^ j — 
, which guarantees the orthogonality of the second 
channel. Channels with F^^ ^ are effectively removed 
from consideration. For example, on the square lattice at 
|;-filling of the upper level, in the absence of the Hund's 
term, $^ would have the highest Tc2 while including the 
Hund's term the calculated Tc2 for this channel turns 
out to be negative. That is no transition would occur if 
channels with F^^ = did not exist. 

To solve for Tc2, it is helpful to isolate the interaction 
contribution to the Green's functions for the delocalized 
level 



(11) 



where E = x x = |Arp<i>^'"(ti^ - \2)^^^\ describes 

the interactions with the localized level. Writing the free 
energy contribution as F\^^\2 = —TTr\n{G~^) + jf ^ 
we find, 



8AraAr ' 



H^=Q — 



N 



2 



EE 

k a=l 



TY,Tr[^i^G",G'M] (12) 

n.k 

fjEka) ~ /(A2) t ^tr^r , 2 

_ VL'^k '^kVka + JT^^ 



where Eka is the energy of the a-th band of the free 
conduction electron problem and we have projected ^j-^fc 
onto the eigenvectors of Hkin', Ai, A2 are set by saddle- 
point evaluations to fix the filling. 

On a square lattice, C4 symmetry means that the pos- 
sible overlap phases with the nearest four neighbors are 
(j)^ = A = {0,..,3}. Hopping occurs in 

the first channel (F = 1) coinciding with the totally sym- 
metric Wannier state A = 0. Time reversal symmetry 
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FIG. 2: Symmetry selects a unique symmetry channel for 
the second-stage Kondo effect in the (a) square lattice (b) 
Kagome lattice and (c)Pyrochlore lattice. The ± signs refer 
to the relative phases of the Wannier state at neigboring sites. 



specifies the d-wave like A = 2 as the second channel 
2). The form factors arise as: J2a'f'*^i'^)'4+a = 



(r 

2 >>S4 



'-i+x ' ^t-y • ^i-x ' ^i+y) which under fourier trans- 
form, defines $^: i/'li — X)fc c^e^^^^^ {cos{kx)+cos{ky)) — 



k ^fe 
*4,3 
k 

The effective T 



Similarly, <i>^ = c{ky) — c{k,j.) are found 



— s{kx) ± is{ky) completing the basis. 



is found from 



E 



/(O) - fiE,p) 



where Ea 



E„ 



-2t|6|2ci>« 



^ = 



(13) 



-|- Ai — fi is determined self- 



consistently subject to the conditions A2 



and 



Efcn/(-2t|6p'i>i')+Ai-M)- 2 • 

On the pyrochlore lattice, 4 bands and 6 nearest neigh- 
bors imply 24 independent symmetry channels! The re- 
quirement that the states preserve the equivalence be- 
tween sites reduces this number to 6, given by 



1 
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(1, a, a^, ±1, ±a, ±a^ 



where a is a cube root of unity. Of these, only the state 
(f)^ — ^(1, 1, 1, — 1, — 1, — 1) is orthogonal to the origi- 
nal channel and time-reversal invariant, uniquely speci- 
fying the second channel. This particular Wannier state 
has the property that it changes sign between neigboring 
tetrahedra on the lattice, as illustrated in Fig. 2. The 
four by four matrix <i>^^^ now has the following parts: 



' c{ka) c{kb) c{kc) ' 

c(ka) c{ka-kh) c{ka-kc) 

c(fcb) c{ka—kij) c{k(, — kc) 

y. c(kc) c{ka — kc) c(fcb~^c) 



(14) 
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while, 



(2) 



S(ka) 
-S(fea) 
— S(fcb) — S(fca— fcb) 



S(fc!,) S(/Cc) 
Si^a — ^b) s{ka—kb) 

— s(fcb — fee) / 



(15) 



For this case, the free hopping Hamiltonian has been 
treated previously by Reimers et al.|Q There are two 
flat bands lying above two dispersing bands. Strictly 
speaking, since we have ^ an electron per site involved 
and four sites, this would actually give us an insulator. 
However, since we have dropped the second degenerate 
band, a more reasonable band structure would likely in- 
clude band crossings and incomplete fillings of the levels 
along the lines of Singh et al. |l7| 

A large body of experimental work has been done, and 
we now discuss a simple interpretation of some of this 
work. We have four localized half-filled electron bands 
interacting via a Kondo-like hybridization with four itin- 
erant levels, each less than quarter-filled, which leaves us 
with a hole- like Fermi surface below Tc2- In this way, our 
model is consistent with the observed sign change in the 
Hall constant which has been observed at 50 K in single 
crystal samples. A simple phenomenological fit consis- 
tent with the ideas we have introduced for the resistivity 
yields 



P = 




T <2K 
T > lOOK 



(16) 



where pi — 21 fiClcm (impurity scattering in the single 
crystal), = 2.0 ^^''™ (arises from the heavy Fermi 
liquid) as given in ref.4, and po = 0.27mf2cTO (impu- 
rity scattering from the localized spin ^ ) , — 1.1 ^^^"^ 
(phonon contribution above Od) and Ap = 2.7*10 
(arising from the 'porcupine' Fermi sea) have been ex- 
tracted from Urano et al.'s data||^. Since ^ ~ (^)^' 

this yields an effective mass ratio of ^ « 27. An anal- 
ogous fit of the magnetic susceptibility yields 



-3 /AHcm 



xT 



yields 



XhT 



T <2K 



2NA{ligYS(S+l)T 
kB{T+6) 



T > 400A 



C17) 



^ « ^ « 29 where the = i (6* = 66K, g 
= 2.21) contribution has been first subtracted away from 
the high temperature susceptibility measurements ||l|, ^, |^, 
^ leaving a contribution linear in T. The observed low 
temperature metal insulator transition is hardly surpris- 
ing under pressure as broken symmetry Wannier states 
lie close in energy to the chosen (jy^ for the second chan- 
nel. Invariance of the heat capacity to SOT is indicative 
that the magnetic coupling J driving the spin quenching 
is of much higher energy than the Kondo temperature 
T«50K. 

In conclusion, we have presented a physical picture 
within which one can understand the physics of the frus- 
trated heavy fermion LiV204. A simple t-J model is 



capable of describing both the high temperature mixed 
valent phase, and the low temperature coupling between 
a Fermi sea and spin liquid. As no time reversal sym- 
metry breaking has been observed, we have been able to 
uniquely identify the symmetry of the low temperature 
quasiparticles. Our picture gives a simple qualitative in- 
terpretation of the unusual resistivity, the change in sign 
of the Hall resistance observed and should be used to pre- 
dict the result of de Haas- van Alphen experiments yet to 
be performed. This work should serve as a starting point 
for a DMFT or band structure calculation which might 
describe the detailed energetics involved. In particular it 
would be interesting to follow the compression of the lat- 
tice and see how the second channel symmetry depends 
on nearest neighbor distance. 
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